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ABSTRACT 



Q ■ We have examined the occurrence of Extremely Red Objects (EROs) in the fields 

^ , of 13 luminous quasars (11 radio-loud and two radio-quiet) at 1.8 < z < 3.0. The 

average surface density of Kg < 19 mag EROs is two-three times higher than in large, 
random-field surveys, and the excess is significant at the ~ 3a" level even after taking into 
account that the ERO distribution is highly inhomogeneous. This is the first systematic 
^ I investigation of the surface density of EROs in the fields of radio-loud quasars above 

■ z ~ 2, and shows that a large number of the fields contain clumps of EROs, similar 

to what is seen only in the densest areas in random-field surveys. The high surface 
densities and angular distribution of EROs suggest that the excess originates in high- 
redshift galaxy concentrations, possibly young clusters of galaxies. The fainter EROs at 
Kg > 19 mag show some evidence of being more clustered in the immediate 20 arcsec 
region around the quasars, suggesting an association with the quasars. Comparing with 
predictions from spectral synthesis models, we find that if the Kg ~ 19 mag ERO excess 
is associated with the quasars at z ^ 2, their magnitudes are typical of > L* passively 
evolving galaxies formed at z ~ 3.5 {^m = 0.3, Q\ = 0.7, and Hq = 70 kms~^ Mpc~^). 
Another interpretation of our results is that the excess originates in concentrations of 
galaxies at z ~ 1 lying along the line of sight to the quasars. If this is the case, the EROs 
may be tracing massive structures responsible for a magnification bias of the quasars. 

Subject headings: galaxies:high-redshift — galaxies: clusters — galaxies: quasars — infrared:galaxies 
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1. Introduction 

The age and space density of the oldest and 
most massive galaxies at each epoch provide 
the strongest constraints on models of galaxy 
formation. This is one of main reasons that 
there has been so much focus on Extremely 
Red Objects (EROs). It is the high-redshift 
systems that are most important for testing 
models, and the ERO selection criterion is de- 
signed to find old, passively evolving galaxies 
at z > 1 . Because the 4000 A break in evolved 
galaxies shifts from optical to near-infrared at 
z >1, EROs are defined in terms of optical to 
near-infrared color, with two of the most used 
definitions being R — K>5.0 and I — K>A.O 
mag. 

The co-moving volume density, or even the 
number density, of massive, evolved galax- 
ies at z > 1 can place strong constraints on 
models of galaxy formation. Whereas mono- 
lithic collapse models (Eggen, Lynden-Bell, & 
Sandage 1962; Larson 1974; Tinsley & Gunn 
1976; Loeb & Peebles 2003) predict pure lu- 
minosity evolution of early- type galaxies from 
the period of formation until the present-day, 
hierarchical formation models (White &: Rees 
1978; Baugh, Cole, & Prenk 1996; Kauffmann 
1996; Baugh et al. 1998; Somerville, Primack, 
&; Faber 2001) predict significant density evo- 
lution from z ~ 1 to z ~ 0. Testing the mod- 
els has proven to be difficult, and the results 
are not yet conclusive. Some surveys have 
found a deficit of early-type galaxies at z ~ 1 
(Barger et al. 1999; Smith et al. 2002; Roche 
et al. 2002) whereas others have not (Daddi, 
Cimatti, k Renzini 2000b; Cimatti et al. 2002; 
Im et al. 2002). Part of the disagreement may 
be caused by the strong clustering of EROs on 
the sky (Daddi et al. 2000a; Firth et al. 2002; 
Roche et al. 2002), and large fields must be 
survoyod in order to overcome this problem. 

^JPL Postdoctoral Researcher 



The ERO definition also selects dusty star- 
burst galaxies whose spectral energy distribu- 
tion can mimic those of old, passively evolv- 
ing galaxies. It is still unclear how large a 
fraction of the ERO population consists of 
dusty starburst galaxies. Several methods 
are being used to separate the two galaxy 
types, including photometric (Mannucci et al. 
2002; Small et al. 2002; Miyazaki et al. 2002), 
spectroscopic (Cimatti et al. 2002; Saracco 
et al. 2003), and those based on morphology 
(Moriondo, Cimatti, k Daddi 2002; Stiavelh 
k Treu 2000; Yan k Thompson 2003). The 
results so far seem to suggest that roughly 
two-thirds of the EROs at i^T < 19 mag are 
old, passively evolving galaxies. However, it 
is also clear that the situation is more com- 
plex than just a simple distinction into old, 
evolved galaxies and dusty starbursts since 
systems with old populations and residual star 
formation also appear in ERO samples (Yan 
k Thompson 2003). 

In this paper, we are concerned with EROs 
in the fields of high-redshift quasars. Some 
early studies serendipitously found EROs 
in the fields of z > 2 radio galaxies and 
quasars (McCarthy, Pcrsson, k West 1992; 
Hu k Ridgway 1994; Dey, Spinrad, k Dick- 
inson 1995). Although some of these EROs 
were later shown to lie in the foreground of 
the radio galaxies (Graham k Dey 1996), 
probing the fields around high-redshift radio- 
loud AGN (active galactic nuclei) has been 
a promising method of finding high-redshift 
galaxies. Luminous quasars and radio galax- 
ies are often associated with high-density 
galaxy environments (Ellingson, Yee, k Green 
1991; Hill k Lilly 1991; Hall k Green 1998; 
Wold et al. 2000; Best 2000; Venemans et al. 
2002; Best et al. 2003), and the rich envi- 
ronments are often interpreted as clusters or 
groups of galaxies hosting the AGN. By com- 
bining a targeted search with a color criterion 
like that of EROs, high-redshift early-type 
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galaxies can be effectively filtered out from 
the numerous foreground galaxies. This be- 
comes very important at high redshifts, since 
even for a rich cluster at z 1, most of the 
galaxies observed using a single broad-band 
filter are foreground or background galax- 
ies (Dickinson 1997). A few surveys have 
found that EROs seem to be more common 
in high-redshift AGN fields than in the gen- 
eral field (Chapman, McCarthy, & Persson 
2000; Cimatti et al. 2000; Hall et al. 2001; 
Sanchez & Gonzalez-Serrano 2002) , and there 
are individual AGN fields which show striking 
overdensities of EROs (Thompson, Aftreth, 
& Soifer 2000; Liu et al. 2000; Haines et al. 
2001). 

Even though galaxy excess in the fields of 
AGN are often being interpreted as groups or 
clusters physically associated with the AGN 
(spectroscopically confirmed in a few cases 
(Dickinson 1997; Deltorn et al. 1997)), gravi- 
tational lensing is also known to cause corre- 
lations (and anti-correlations) between high- 
redshift quasars and optically bright fore- 
ground galaxies (Webster et al. 1988; Hammer 
& LeFevre 1990; Benitez, Martinez-Gonzalez, 
& Gonzalez-Serrano 1995; Benitez, Martinez- 
Gonzalez, &; Martin-Mirones 1997; Norman & 
Impey 2000; Myers et al. 2003). It is possible 
that a similar situation could arise between 
high-redshift quasars (z > 2) and EROs, be- 
cause the EROs are likely to be tracers of mas- 
sive structures which may boost the fluxes of 
distant background quasars by gravitational 
lensing. 

The aim with our present study is to test 
whether EROs are more common in the fields 
of z > 2 quasars. To our knowledge, this is 
the first systematic survey to probe the sur- 
face density of EROs in the fields of quasars at 
z > 2, and preliminary results were reported 
by Wold et al. (2003). The quasars in our sam- 
ple are selected from the catalogs of Barthel et 
al. (1988) and Hewitt k Burbidge (1993), and 



form a subset of a larger sample of 40 designed 
to study host galaxies (Armus et al. in prepa- 
ration). The 40 sample quasars are nearly 
evenly divided between radio-loud and radio- 
quiet systems with comparable V magnitudes 
and redshifts. The subsample of 13 analyzed 
here were selected randomly from the sample 
of 40 and followed up with wide-field near- 
infrared and optical imaging. The sample is 
listed in Table 1. In this table, we have also 
listed the radio loudness, defined in terms of 
the ratio of the flux densities at rest-frame 
5 GHz and 2500 A, S^ghz/ S2500 (Sramek & 
Weedman 1980). According to the definition 
by Kellermann et al. (1989), a quasar is radio- 
loud if S^GHz/ S2500 > 10) and thus two of the 
13 quasars arc classified as radio-quiet (one of 
these might be radio-intermediate instead of 
radio-quiet) whereas the rest are radio-loud. 
In calculating the radio loudness, we used the 
i^s-band fluxes of the quasars as found from 
our data, and the radio fluxes at 5 GHz listed 
in the catalog of Barthel et al. (1988) unless 
otherwise is noted in the table. The Kg-hand 
fluxes were converted to rest-frame 2500 A by 
assuming a power-law spectrum, oc z^~", 
where a was found from the median quasar 
energy distribution by Elvis et al. (1994). The 
conversion to rest-frame 5 GHz radio flux was 
done by assuming the same power-law form of 
the spectrum, but using spectral indices listed 
in the catalog by Barthel et al. (1988), unless 
otherwise is noted in the table. 

The structure of the paper is as follows. In 
the next section, we describe the observations 
and data reduction. Section 3 deals with the 
construction of object catalogs and addresses 
their incompleteness. In Section 4 we present 
the results of the analysis. The surface den- 
sity of EROs is calculated and the ERO ex- 
cess evaluated, carefully taking into account 
their strong clustering. We also investigate 
the i^'^-band number counts and the radial 
distribution of EROs around the quasars. In 
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Table 1 
The quasar sample. 



Quasar 


SsGHz 
52500 


z 


R 


Ks 


(1) 


(2) 


(3) 


(4) 


(5) 


0741+169 


3880^ 


1.894 


7500 


3280 


0758+120 


1690 


2.660 


3600 


2360 


0830+115 


2050 


2.974 


4550 


2720 


0927+217 


130 


1.830 


3600 


3280 


0941+261 


5750 


2.910 


2400 


2240 


1056+015 


<14^ 


2.650 


11100 


3260 


1354+258 


1080 


2.032 


2500 


6800 


1456+092 


2590 


1.991 


3900 


3200 


1629+680 


4070 


2.475 


6100 


3780 


1630+374 


<6'' 


2.037 


3600 


2920 


1702+298 


18130 


1.927 


4000 


3340 


2150+053 


3330 


1.979 


2000 


3040 


2338+042 


4915 


2.594 


3600 


2880 



Note. — (1) lAU designation. (2) Radio loudness defined as the ratio of the flux densities at 
rest-frame 5 GHz and 2500 A (Sramek k Weedman 1980). (3) Redshift. (4)-(5) Total exposure 
time in seconds in the R- and i^^-filters. 

'^Radio Atix density and spectral index from the Green Bank radio survey by Becker, White, & 
Edwards (1991). 

^Upper limit based on a 3a detection in the NRAO-VLA Sky Survey (Condon et al. 1998) and 
an assumed spectral index of a = 0.7. 
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Section 5, we discuss two different interpreta- 
tions of the results, and the conclusions are 
drawn in Section 6. 

All magnitudes are given in the Vega sys- 
tem, and our assumed cosmology has = 
0.3, = 0.7, and Ho = 70 kms"^ Mpc^^ 

2. Observations and data reduction 

The data consist of wide-field R- and Kg- 
band images taken with the Palomar 200-inch 
telescope using the prime-focus instruments 
COSMIC^ (Kells ct al. 1998) and PFIRCAM^ 
(Jarrett et al. 1994a). The COSMIC instru- 
ment is equipped with a 2048x2048 Tektronix 
CCD with a pixel scale of 0.2856 arcsec, and 
the PFIRCAM employs a 256x256 HgCdTe 
array with a pixel scale of 0.494 arcscc. The 
resulting fields of view in R and Ks are there- 
fore 9.7x9.7 and 2.1x2.1 arcmin^, respec- 
tively. 

Images were obtained during several ob- 
serving runs from September 1995 to July 
2001. We used a 9- or 27-point dither pat- 
tern for the Kg images, and on-chip integra- 
tion times of 4 sec. The number of co-adds 
was set so that the total integration time at 
each dither point was 20-30 sec, and the off- 
sets between individual frames were typically 
10-15 arcsec. The i?-band exposures were di- 
vided into integrations of 300-600 sec each, 
with an offset of 10-20 arcsec between each 
exposure. The total integration time in R and 
Kg for each field is listed in Table 1. 

Flat-field images in Kg were made by tak- 
ing the median of dark-subtracted quasar im- 
ages. Typically, a flat-field was made for each 
dither sequence, and this was found to be suf- 
ficient in most cases. After flat-fielding, the 
images were sky subtracted, corrected for bad 
pixels and registered using the DIMSUM pack- 

^ Carnegie Observatories Spectrograph and Multiple 
Imaging Camera 
^Prime Focus Infrared Camera 



age (Stanford, Eisenhardt 8z Dickinson 1995) 
within IRAF.^ 

The sky subtraction in DiMSUM is per- 
formed by subtracting a running sky frame 
made by taking the median of nine neighbor- 
ing images. After the sky subtraction, masks 
are created for objects and cosmic rays and 
the masks are utilized to make a first pass 
combined image. The first pass image is then 
used to make a new and improved mask by in- 
cluding fainter objects. The new object mask 
is used in the final pass where the sky sub- 
traction is repeated and the images aligned 
and combined. At this final step, we block 
replicated the images by a factor of two so 
that sub-pixel shifts were used for the align- 
ing. The resulting images therefore have a 
pixel scale of 0.247 arcsec. 

The i?-band images were reduced in IRAF 
in a standard manner. The bias level was sub- 
tracted either by using a bias frame or by us- 
ing a 3rd order Legendre to fit the overscan 
region which was then subtracted from each 
row. A flat-field was formed for each quasar 
field by taking the median of the quasar im- 
ages, or of the quasar images in combination 
with exposures of other fields taken near in 
time. Typically 15-30 frames were used for 
this, and all emission ?a 2a above the average 
background was masked out. This method re- 
jected cosmic rays, stars and galaxies visible 
in the quasar images and yielded a high qual- 
ity flat-field frame. The flat-fielded images 
were then aligned using 20-40 stars in each 
frame, and averaged using imcombine with a 
high threshold set to remove bad pixels that 
had been flagged by giving them high, dis- 
crepant values. The combined R images were 
aligned with, and put to the pixel scale of, the 

*IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association 
of Universities for Research in Astronomy, Inc., un- 
der cooperative agreement with the National Science 
Foundation. 
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Kg images using the IRAF tasks geomap and 
geotran. 

The flux cahbration of the i?-band images 
was accomphshed using standard stars from 
Landolt (1992), and the Kg-daXa were cah- 
brated using standard stars from Pcrsson ct 
al. (1998). The conditions were photometric, 
and the transparency good throughout most 
of the runs. The final photometric status was 
checked by the amount of variation within 
the individual exposures of the same source 
and by following standard stars throughout 
the night. For fields that were imaged during 
non-photometric conditions, calibration im- 
ages were obtained later during photometric 
conditions to allow for an accurate flux cali- 
bration. During the photometric nights, the 
zeropoint varied by less than 0.1 mag. For a 
subset of the images, we checked our absolute 
photometry against the Two Micron All Sky 
Survey (2MASS) point source catalog (Skrut- 
skie et al. 1997; Cutri et al.^). A very good 
agreement, typically < 0.1 mag difference, 
was found between our Ks-magnitudes and 
the 2MASS if^-magnitudes for bright stars of 
Kg < 15 mag. 

The full width at half maximum (FWHM) 
of the seeing in the X^-images is 0.8-1.0 arc- 
sec, and in the i?-images, typically 1.0-1.5 
arcsec. The raw number counts of galaxies 
with detection significance > 3cr turns over at 
Ks ~ 20 mag, although this varies somewhat 
from field to field due to the different expo- 
sure times. In the i?-band the corresponding 
turn-over in the number counts is i? ~ 24.5 
mag. The 3a detection limits in the images 
are Kg ~ 21.5 and R ^ 26 mag. We address 
the completeness of the images in the next 
section. 

® http : // www.ipac. calt ech . edu/ 2mass /releases /second/ 
doc/ explsup .html 



3. Construction of object catalogs 

3.1. Source detection 

For object detection and photometry we 
used the SEXTRACTOR software package (Bertin 
& Arnouts 1996). In order to have the same 
FWHM of the seeing in both R and Kg we 
convolved the Kg-hand images with a Gaus- 
sian function using the IRAF task gauss prior 
to object detection. For each field, object 
catalogs were constructed based on the detec- 
tions in the i^^-band image. During the detec- 
tion process, we used a Gaussian convolution 
filter with a FWHM matched to the seeing 
and a detection threshold of 1.50". Since the 
images were block replicated, we set the mini- 
mum number of pixels above the threshold to 
20 in order to avoid spurious detections and 
cosmic rays over an area of less than 5 pixels 
in the original images. The result of the detec- 
tion process was inspected visually in order to 
ensure that no obvious objects were missed, 
and that no false detections were entered into 
the catalogs. Saturated objects and objects 
lying close to the image boundaries were re- 
jected from the catalogs. 

Because of the dithering technique used in 
the Kg-hand, the exposure time falls off to- 
ward the edges in the final combined images. 
Since the PFIRCAM field of view is relatively 
small, we included lower exposure regions, 
but rejected the regions where the exposure 
time is less than 25 % of the total, i.e. where 
the signal-to-noise drops by more than 50 %. 
This causes the rms background noise to in- 
crease toward the edges, resulting in a num- 
ber of spurious detections. In order to avoid 
such spurious detections, we used the expo- 
sure maps produced by DIMSUM as weights 
in SEXTRACTOR. This allows for an adjust- 
ment of the detection threshold according to 
the noise level. 

After the i^^s-catalogs were made, matched 
i?-band catalogs were created by running SEX- 
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14 16 18 20 22 

Ks [mag] 

Fig. 1. — The circles show objects detected 
at a significance > 5a in the i^Ts-filter. The 
triangles are objects with a detection signifi- 
cance < 3a in the i?-filter where a lower limit 
to the R — Kg color has been calculated as- 
suming a 3a detection in R. The crosses are 
objects classified as stars. The horizontal line 
indicates our ERO definition, and the diago- 
nal line the R = 26 mag limit. The scatter 
of objects above the diagonal line arises from 
the varying depths of the i?-band images. 



TRACTOR in double image mode. This pro- 
duces a catalog of fluxes measured in the 
R-hand image, but using the i^^-band im- 
age for detection. For total magnitudes, we 
used the sextractor best magnitudes which 
are evaluated using adaptive apertures, or us- 
ing isophotal corrections if a neighbor is sus- 
pected of biasing the adaptive aperture mag- 
nitude (Bertin Sz Arnouts 1996). For deter- 
mining colors, we used aperture magnitudes, 
and chose an aperture diameter of two times 
the FWHM of the seeing, corresponding to 
typically 9-10 pixels, or 2.2-2.5 arcsec. 

For each ETg-selected object, we calculated 
the detection significance using the uncer- 
tainty in the flux as measured by sextrac- 
tor which takes into account the Poissonian 
nature of the counts and the standard devia- 
tion in the background counts. In addition to 
this, we added in quadrature the error arising 
from the fluctuating background level across 
the image (Best 2000). In order for a Kg- 
selected object to be included in the catalog 
we required the detection significance within 
the aperture to be > 5a. If the flux measured 
within the corresponding aperture in the R- 
band was significant at the > 3a level, we 
evaluated the R — Kg color, and if < 3a, a 
lower limit to the R — Kg color was calculated 
assuming a 3a detection in R. 

In order to separate stars and galaxies, we 
used the star-galaxy classifier in sextractor 
which assigns a probability between and 1 
to an object according to how likely it is to 
be a star. The star-galaxy classifier works 
well, but becomes less reliable at faint mag- 
nitudes. We therefore performed star-galaxy 
separation only at Kg < 17 mag, rejecting ob- 
jects with a star-galaxy class of > 0.85 from 
the catalogs. The star-galaxy separation was 
based on the detections in the Kg images. 
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3.2. Completeness 

To be able to address the completeness 
of the ii's-band catalogs, wc added artifi- 
cial galaxies to the images using the IRAF 
task mkobjects. For a given input magni- 
tude, we added galaxies with both deVau- 
coulcurs and exponential disk profiles, in to- 
tal ~ 1500 galaxies distributed over the 13 
fields. The scale- lengths of the artificial galax- 
ies were randomly selected from the FWHM 
of the objects in the catalogs, and the images 
convolved with the same Gaussian function as 
was used to put the i^^-band images to the 
seeing of the i?-band images. After this, the 
images were processed through SEXTRACTOR 
using the same detection criteria as described 
above. 

The simulations show that for input ob- 
jects of Kg = 18 mag, 96 % of the galaxies 
are recovered, where the missing 4 % is due 
to blending with other objects. At Kg = 19 
mag the completeness drops to 90 %, and at 
Kg = 20 mag the catalogs are 40 % com- 
plete. The numbers were evaluated assuming 
a 50:50 mix of deVaucouleurs and exponential 
disk profiles. 

4. Results 

4.1. The surface density of EROs 

We select EROs in the quasar fields using 
the criterion R — Kg > 5.0 mag. This includes 
old, passively evolving galaxies at z > 0.9, but 
also dusty starbursting galaxies. With only 
two filters, and at the resolution of our im- 
ages, it is not possible to distinguish between 
these two galaxy types. There will also be 
some contamination by faint, red stars in our 
ERO sample as discussed below. We have cho- 
sen the R — Kg > 5.0 mag selection criterion 
because of two reasons. Firstly, because se- 
lecting galaxies with red colors enables us to 
eliminate a large part of the numerous fore- 
ground galaxies at z < 1 while at the same 



time producing a big enough sample to per- 
form meaningful statistics on. Secondly, be- 
cause several random-field surveys employing 
the R — Kg > 5.0 mag criterion exist with 
which we can compare the ERO surface den- 
sity. Our areal coverage in R and Kg is limited 
by the small field of view of the PFIRCAM 
instrument, so the field counts cannot be es- 
timated from the quasar images themselves. 

The color-magnitude diagram of all the Kg- 
selcctcd objects is shown in Fig. 1. Above our 
R — Kg > 5.0 mag selection there are 253 ob- 
jects, and 105 of these have a detection sig- 
nificance in R of three or less, i.e. the lower 
limit on their R — Kg color is > 5.0 mag. The 
SEXTRACTOR star-galaxy classifier finds that 
two of the EROs at Kg < 17 mag are stars, 
and these are excluded from the ERO sample, 
resulting in a total of 251 galaxies (including 
the < 3(7 detections in R having a lower limit 
oiR-Kg> 5.0 mag). 

In Fig. 2, we have plotted the average 
cumulative surface density of EROs in the 
quasar fields together with ERO surface den- 
sities found in the random-field surveys of 
Daddi et al. (2000a), Scodeggio & Silva 
(2000), Cimatti et al. (2002), Roche et al. 
(2002), and Miyazaki et al. (2002). Note that 
we have not compared with field counts de- 
rived from other colors (e.g. I — K or I — H) 
since they may bias toward later-type galax- 
ies with prolonged star formation (McCarthy 
et al. 2001; Yan & Thompson 2003). Note 
also that Miyazaki et al. define EROs as hav- 
ing R — Kg{AB) > 3.35 mag, which, ac- 
cording to Miyazaki et al., corresponds to 
R - Kg{Vega) > 4.95 mag. 

The literature surveys cover wide, random 
areas and should therefore be representative 
of the ERO counts in the general field. The 
field counts agree well with each other ex- 
cept for the counts in the Chandra Deep Field 
South by Scodeggio & Silva (2000), which are 
consistently lower than the other. This is 
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most likely caused by the inhomogcncity of 
the ERO distribution coupled with the rela- 
tively small (and contiguous) area surveyed 
by Scodeggio & Silva. Fig. 2 demonstrates 
that the ERO counts in our quasar fields arc 
systematically higher than those in the field. 
At Kg < 19 mag, where our counts are com- 
plete, the average surface density of EROs is 
two-three times larger than the average of the 
field values (excluding those in the Chandra 
Deep Field South). According to the average 
of the random-field surveys, as indicated by 
the dotted line in Fig. 2, we expect 34 it 14 
EROs at Kg < 19 mag over an area equal 
to the total area of our survey. (The error 
in this number takes into account the strong 
clustering of EROs as explained below.) The 
observed number of EROs in the quasar fields 
at < 19 mag is 79, hence there is a signifi- 
cant excess above the field counts. At Kg > 20 
mag we still see a clear excess in the quasar 
fields, but our counts become strongly affected 
by incompleteness. 

The number of EROs per quasar field and 
the average cumulative surface density as a 
function of ii's-magnitude is listed in Table 2. 
Some fields are seen to be very rich in EROs, 
such as 1456+092 and 2150+053, with surface 
densities four-five times larger than the field 
values. These fields have ERO excesses in ev- 
ery magnitude bin. Other fields have surface 
densities consistent with the field counts, like 
0741+169 and 1629+680. The majority of the 
fields, however, have surface densities in be- 
tween these two extremes. This is illustrated 
in Fig. 3, where a histogram of the ERO sur- 
face densities at Kg < 19 mag is shown. If 
the quasar fields were random with respect 
to the ERO distribution, we would expect 
a normal distribution centered on ~ 0.546 
arcmin"^, which is the average surface density 
in the general field. This is not seen, instead 
there is a tail toward high surface densities, 
with as many as four fields (30 per cent of 



the sample) having a surface density larger 
than three times the field value, and seven 
fields where the surface density is twice the 
field value. The combination of the very rich 
and the moderately rich fields gives an overall 
surface density well above the average for the 
random-field surveys. 

4.2. Contamination by red stars 

Of the 13 quasars fields, 10 lie at high galac- 
tic latitudes b > 35°, so the contamination by 
stars is likely to be small. Using a color-cut 
of R — Kg > 5.0 mag will eliminate most field 
stars because of their bluer color, but L- and 
M-type dwarfs will be red enough to make it 
into our ERO sample. Based on R — K' and 
J-K' colors, Mannucci et al. (2002) find that 
9 % of their K' < 20 mag ERO sample (at 
b = —23°) is likely to consist of stars. 

In order to estimate the contamination by 
stars, we used the stellar distribution model 
developed by Jarrett (1992) (see also Jarrett, 
Dickman, k Hcrbst (1994)). The model is 
an extension of the Bahcall & Soneira (1980) 
optical star count model to the near-infrared, 
and its performance with respect to 2MASS 
near-infrared star counts is demonstrated by 
Cambresy et al. (2002). For every quasar po- 
sition, we used the model to predict the ex- 
pected number of stars having R — Kg > 5.0 
mag. The predicted average surface density 
of stars for the whole ERO sample is listed 
on the bottom line of Table 2. These num- 
bers show that 40-50 % of the brighter EROs 
at Kg < 18 mag might be stars, but that 
at fainter magnitudes the contamination de- 
creases to 15-20 %. Note that this is an over- 
estimate because the ERO counts suffer from 
incompleteness at fainter magnitudes. The 9 
% contamination estimated by Mannucci et 
al. (2002) based on two colors therefore seems 
to agree well with the model predictions. We 
thus expect roughly 10 % of our EROs to 
be red stars. Since we cannot use the small 
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Fig. 2. — The cumulative surface density of EROs in the quasar fields compared with literature 
surveys in R — K or R — Kg covering wide, random fields (Roche et al. 2002 (R02); Cimatti et al. 
2002 (C02); Miyazaki et al. 2002 (M02); Scodeggio k Silva 2000 (SSOO); Daddi et al. 2000a (DOO)). 
The average of the literature values (excluding the SSOO counts) is indicated by the dotted line. 
Error bars are Poissonian. 
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Table 2 

The cumulative number of EROs in the quasar fields as a function of total K, 

magnitude. 



villas ax 


Area 


^ icS.U 


^ ±0.0 


<' 1 Q n 
^ ±y.u 


^ ±y.o 


^ on n 


<' on ^ 




arcmin^ 


mag 


mag 


mag 


mag 


mag 


mag 


0741+169 


5.22 





1 


4 


7 


9 


10 


0758+120 


5.23 


2 


5 


6 


8 


12 


16 


0830+115 


5.28 


1 


3 


7 


10 


14 


18 


0927+217 


4.19 


1 


4 


8 


12 


18 


22 


0941+261 


4.16 





2 


5 


8 


14 


17 


1056+015 


5.13 





2 


4 


6 


8 


16 


1354+258 


5.07 





3 


5 


10 


14 


17 


1456+092 


4.34 


2 


4 


10 


12 


15 


19 


1629+680 


5.09 


2 


3 


4 


11 


13 


21 


1630+374 


4.30 


1 


2 


2 


5 


5 


5 


1702+298 


5.13 


1 


1 


4 


9 


14 


15 


2150+053 


4.21 


2 


7 


10 


13 


16 


21 


2338+042 


4.36 


2 


6 


10 


16 


21 


27 


N 

qso 


61.7 


14 


43 


79 


127 


173 


224 


Uqso (arcmin"^) 




0.23+0.06 


0.70+0.11 


1.28+0.14 


2.06+0.18 


2.80+0.21 


3.63+0.: 


nstars (arcmin~2) 




0.11+0.04 


0.17+0.05 


0.25+0.06 


0.44+0.08 


0.47+0.09 


0.59+0.: 



Note. — The mimbcrs in this table are based on raw counts. The bottom three Unes show the 
total number of EROs in the 13 fields, Ngso, the resulting average surface density, Ugso, with its 
associated Poissonian error, and the expected surface density of stars, Ugtars, with its associated 
Poissonian error. According to random-field surveys, the expected number of < 19 mag EROs 
over 5 arcmin^, the typical size of one quasar field, is two-three. 
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quasar fields for estimating the field counts of 
EROs, we compare the counts in the quasar 
fields with number densities of EROs from 
random-field surveys found in the literature. 
Since the random-field STirvcys arc also con- 
taminated by stars, a contamination of 10-20 
% does not alter our conclusions. 

The three lowest galactic latitude fields 
in our sample are 0741-1-169, 0758-1-120 and 
0830-^115 at 6 = 20 - 27°. For these fields 
the model predicts 2.9, 2.6, and 1.6 stars at 
Kg < 19 mag, respectively. Note also that 
these are not among the richest fields in our 
sample. 

4.3. Significance of the ERO excess 

In this section we evaluate the significance 
of the ERO excess, taking into account the 
large variance in the ERO field counts caused 
by their strong clustering. The significance of 
the excess is given by 



i^qso - Afield) 



a 



(1) 



(Yee & Lopez-Cruz 1999), where Nqgo is the 
total number of EROs in the quasar fields 
brighter than a given i^^-magnitude and 
Afield is the number of EROs in the random 
field scaled to the total area of the quasar 
fields, 61.7 arcmin^. The error in the field 
counts expected over an area of 61.7 arcmin^ 
is denoted a. Ideally Nfi^i^ should have been 
estimated from the edges of the quasar fields, 
but our fields are too small for this. Instead 
we use the average of the surface densities 
found by Daddi et al. (2000a), Roche et al. 
(2002), and Miyazaki et al. (2002). 

Because ERO positions are correlated on 
the sky, the variance is given by 



(t'^ = N field (1 + NfieldA^C) 



(2) 



(Daddi et al. 2000a), where A^^ is the ampli- 
tude of the ERO angular two-point correla- 
tion function. The integral constraint which 



arises from the fact that the survey area is 
limited is given by A^C, and Daddi et al. find 
that C can be approximated by 58 x Area~^'^. 
We use this approximation and find that C = 
11.15 for the 13 quasar fields. The approxima- 
tion is valid for a contiguous area, but serves 
our purpose well even if the 13 fields are non- 
contiguous. In any case, C = 11.15 is likely to 
be an overestimate since our fields are widely 
separated and not correlated in any way. We 
use the amplitudes, Ai_j, found by Daddi et 
al. and Roche et al. for Kg < 19 mag and by 
Roche et al. for Kg < 19.5 mag, and calculate 
a. 

The result of the calculations for the entire 
quasar sample on average is shown in Table 3, 
where it can be seen that the excess at Kg < 
18, 18.5 and 19 mag is significant at the So- 
level. Because of incompleteness, the excess 
is not significant at fainter i^Ts-levels, but if a 
correction is made for the incompleteness, the 
excess appears to be significant also at fainter 
Kg magnitudes, as shown in the last column 
of Table 3. 

4.4. Radial distribution 

If the EROs are associated with the quasars 
either through a surrounding cluster or a mag- 
nification bias we might expect to see a trend 
with projected quasar separation. In order 
to examine the radial distribution of EROs, 
we evaluated the average surface density in 
three annuli centered on the quasars, for three 
different magnitude intervals. The result is 
shown in Fig. 4. The brighter EROs at 
17.5 < Kg < 19.5 mag are seen to be uni- 
formly distributed over the whole field, but 
the fainter ones at Kg > 19.5 mag seem to 
cluster more in the immediate 20 arcsec re- 
gion around the quasars. Note that this is 
not caused by the mosaicing technique in the 
Kg filter. The exposure times across the inner 
two annuli are uniform, with the third annulus 
being on the average only slightly less well ex- 
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Table 3 

Significance of the average ERO excess in the quasar fields. 



(1) 


Nqso 

(2) 


N field ± O- 

(3) 


Nexc 

(4) 


^exc 

(5) 


-.corr 
^ exc 

(6) 


< 18.0 


14 


4.9±3.4 


9.1 


2.7 




< 18.5 


43 


17.2±9.5 


25.8 


2.7 


3.0 


< 19.0 


79 


33.7±14.3 


45.3 


3.2 


4.0 


< 19.5 


127 


72.8±31.1 


54.2 


1.8 


2.5 


< 20.0 


173 


111.9±34.4 


61.1 


1.8 


3.5 


< 20.5 


221 


153.1±1G.1 


70.9 


1.5 





Note. — (1) Ks magnitude limit. (2) Total number of EROs in the 13 quasar fields. (3) The 
expected total number in the 13 fields based on the average of the literature counts, where a is the 
uncertainty taking into account that the ERO distribution is inhomogeneous. (4) Excess above the 
field counts. (5)-(6) Significance of the ERO excess, calculated using Eq. 2. All numbers, except 
those in column (6), were evaluated using raw counts. 
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Fig. 3. — Histogram of raw ERO surface den- 
sities at Kg < 19 mag in the quasar fields. 
The vertical, dashed line indicates the aver- 
age of the random-field surveys, n = 0.546 



arcmm 



posed than the other two. Using completeness 
simulations wc find that the variation in the 
completeness level across the three annuli is 
at most a few per cent, i.e. consistent with 
what we expect from object crowding and 
blending. Therefore, there appears to be a 
real trend toward central concentration of the 
fainter EROs. Cimatti et al. (2000) looked for 
central concentrations of i^s-selected galaxies 
toward z f« 1.5 radio-loud AGN. They did 
not find any significant trends, but noted that 
the regions closest to the AGNs had a higher 
number density of galaxies than in the general 
field. 



4.5. The surface density of R — K^ > 6.0 
mag EROs 

For completeness, we also evaluated the 
surface density oi R — > 6.0 mag EROs. 
Using the redder color criterion gives a much 
smaller sample, 14 EROs at Ks < 19 mag 
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Fig. 4. — Mean radial distribution of EROs 
in the 13 quasar fields based on raw counts. 
The inner bin is centered on the quasars. The 
error bars show the standard deviation in the 
mean. 



and 31 EROs at Kg < 20 mag. The result- 
ing surface densities are thus 0.23 it 0.06 and 
0.50 ± 0.09 arcmin-2 at Kg < 19 and < 20 
mag, respectively (Poissonian errors). The 
surface density of EROs at Kg < 19 mag is 
therefore 3-5 times higher than that found in 
the field surveys by Thompson et al. (1999) 
and Daddi et al. (2000a) (0.039 ± 0.016 and 
0.07 lb 0.01 arcmin"^, respectively). We also 
note that the surface density of R — Kg > 6.0 
EROs in our quasar fields agrees with the 
ERO surface densities found in radio-loud 
AGN fields at z = 1 - 1.5 by Hall & Green 
(1998) and Cimatti et al. (2000). 

4.6. Ks-band counts 

Whereas there is a clear excess of EROs 
in the quasar fields, the overall if^-counts 
show only a slight excess at fainter magni- 
tudes, Kg > 18.5 mag. This is seen in Fig. 5 
where we have plotted the average i^s-band 
number counts together with counts from the 
surveys by Soifcr ct al. (1994), Moustakas 
et al. (1997), and Maihara et al. (2001), as 
well as the literature averages compiled by 
Hall, Green, & Cohen (1998) and Best et al. 
(2003). There are large field-to-field varia- 
tions in the random-field i^g-band counts, and 
even though the counts in the quasar fields 
appear to lie above the literature averages at 
fainter Kg magnitudes, they are still consis- 
tent with the counts from the surveys by e.g. 
Soifer et al. (1994), Moustakas et al. (1997), 
and Maihara et al. (2001). We also find that 
the i?-band number counts are consistent with 
the counts in the general field (Metcalfe et al. 
1995; Hogg et al. 1997). Overplotted in Fig. 5 
are also the number counts of EROs in our 
fields, demonstrating, as expected, that non- 
ERO galaxies clearly dominate at all magni- 
tudes. 

Table 4 shows the fraction of EROs with 
respect to the whole if^-selected sample of 
galaxies. For comparison, we have also listed 
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the ERO fractions given by Daddi ct al. 
(2000a), Scodeggio & Silva (2000), and Roche 
et al. (2002), the latter calculated from their 
table 2. In the quasar fields, the EROs make 
up typically 20-25 % of all the ii's-selected 
galaxies, whereas the field surveys find < 10 
%. We note that this behavior is what we 
expect if the EROs are tracing the early-type 
population in high-redshift clusters of galax- 
ies since the fraction of early-type galaxies 
in rich clusters is known to be different from 
the field (Dressier 1980). This also illustrates 
the richness of EROs in the quasar fields, and 
suggests that the EROs arc efficiently trac- 
ing overdensities not readily apparent from 
the Kg number counts alone. Clearly, the 
color-selection is very effective in eliminating 
foreground contamination. 

5. Discussion 

We have demonstrated above that the 
quasar fields have a significant excess of EROs 
above what is expected from the general field 
population of EROs, taking into account their 
strong clustering. It appears that the quasar 
fields contain clusters or dense concentrations 
of red galaxies, and their radial distribution 
suggests that the fainter EROs may be related 
to the quasars. However, another possibility is 
that the EROs are related to structures in the 
foreground of the quasars. We discuss each 
of these possibilities below. Note that our 
sample is dominated by radio-loud quasars, 
and that we restrict the following discussion 
to radio-loud quasars only. 

5.1. 'Clusters' of EROs associated 
with the quasars 

An argument in favor of the ERO excess 
being at the quasar redshifts is that lumi- 
nous radio-loud quasars and radio galaxies 
at z < 1 are known to reside in regions of 
enhanced galaxy density, interpreted as clus- 



ters or groups of galaxies physically associ- 
ated with the AGN (Hall & Green 1998; Best 
2000). Since powerful AGN with supermas- 
sive black holes may trace the densest regions, 
it is natural to assume that this is true also 
at z > 2. For example, Cimatti ct al. (2000) 
find excess EROs in radio-loud AGN fields at 
z ^ 1.5 and Best et al. (2003) find that red 
galaxies cluster around powerful radio sources 
at z 1.5. Our study pushes this out to 
z>2. 

If the ERO excess originates in dense en- 
vironments associated with the quasars, the 
majority of the EROs must lie at z ~ 2 since 
this is the rcdshift of most of the quasars in 
our sample. It is therefore appropriate to dis- 
cuss whether the color and the brightness of 
the EROs are consistent with being evolved 
galaxies at z ^ 2. We have used the pegase 
models of Fioc & Rocca-Volmerange (1997) 
to examine both the apparent magnitude and 
the R — Kg color of passively evolving galax- 
ies at z ^ 2 (assuming no extinction by dust, 
and an initial metallicity of 1 % solar value). 
The apparent Kg magnitude of L* galaxies at 
z = 2 was found by calculating the scale fac- 
tor between the observed flux of the model 
and that of a present-day L* galaxy. The 
present-day L* galaxy was assumed to have 
= -20.3 mag (Loveday et al. 1992) and 
an age of 14 Gyr. 

The standard model in which the stars 
form in a 1 Gyr long burst at ^ > 4 has an 
R — Kg color at z = 2 which is too red for 
most of our EROs. In order for this model 
to produce a color oi R — Kg = 5 — 6 mag 
at 2; = 2, the formation redshift has to be 
z 3.5. The K* magnitude of such a galaxy 
observed at z = 2 is 19.2 mag. Since the star 
formation in this model stops abruptly after 1 
Gyr, this model represents the galaxy which 
reddens most quickly. Therefore, Kg = 19.2 
mag is the brightest magnitude we can expect 
for an L* galaxy at z = 2. If the Kg « 19 
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Fig. 5. — Number counts of i^^s-selected galaxies (circles) and EROs (diamonds) detected at a 
significance greater than 5. We have compared with a selection of deep i^T-band surveys from the 
literature as shown by the other symbols. The dotted and dashed lines show the literature averages 
compiled by Hall, Green, &; Cohen (1998) and Best et al. (2003), respectively. 
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Table 4 

ero fractions among ks-selected galaxies. 



Ks 


This work 


DOO 


SSOO 


R02 


(1) 


(2) 


(3) 


(4) 


(5) 


< 18.0 


0.14±0.04 


0.05 






< 18.5 


0.22±0.04 


0.08 






< 19.0 


0.25±0.03 


0.12 


0.02±0.01 




< 19.5 


0.26±0.03 




0.06±0.02 


0.12±0.02 


< 20.0 


0.24±0.02 




0.06±0.01 


O.MibO.Ol 


< 20.5 


0.22±0.02 









Note. — (1) Magnitude limit. (2) Fraction of EROs evaluated using raw counts. (3)-(5) ERO 
fractions found by Daddi et al. 2000a (DOO), Scodeggio & Silva 2000 (SSOO), and Roche et al. 2002 
(R02). The uncertainties are Poissonian. 



mag ERO excess is associated with passively 
evolving galaxies at 2; a; 2, they must therefore 
be > L* galaxies. Other models which have 
R — Kg = 5 — 6 mag at z = 2 include those 
with exponentially decreasing star formation 
rates. For instance, a galaxy formed at z = 6 
with an exponentially decreasing star forma- 
tion rate of r = 0.5 Gyr has R — Ks = 5 — 6 
mag at z = 2, and = 19.6 - 19.7 mag. If 
the EROs are described by this latter model, 
they must have small amounts of residual star 
formation in them. 

The recent spectroscopic work of van 
Dokkum et al. (2003) has shown that red 
galaxies at 2; > 2 can be very bright. Five 
of their six galaxies were found to lie at 
2.4 < z < 3.5 and to have i^^-magnitudes in 
the range 19.2-19.9 mag. These galaxies were 
selected by the J — Kg > 2.3 mag criterion 
(Franx et al. 2003), whereas our galaxies are 
selected on the basis of having R — Kg > 5.0 
mag. Both criteria are designed to find galax- 
ies with evolved 4000 A breaks, but at dif- 



ferent redshifts. However, these two selec- 
tion techniques will not necessarily sample 
the same type of galaxies since the different 
filters sample different rest-frame wavelength 
ranges. While the results of van Dokkum ct 
al. prove that high-redshift red galaxies can 
be bright, their spectroscopic sample (six of 
eleven J— i^^-selected galaxies) may be biased 
toward systems with ongoing star formation. 

The average ERO fraction in our quasar 
fields is above 20 %, and for the four richest 
fields in our sample the fraction is as high as 
~ 50 %. This is what we expect if the fields 
contain clusters or concentrations of galax- 
ies with an already existing relation between 
morphology and density (Dressier et al. 1997). 
This is also the suggestion from the work by 
Best et al. (2003) who find that red galaxies 
are concentrated toward 2; 1.5 powerful ra- 
dio sources. 

Whereas the magnitudes and colors of the 
EROs are consistent with models of passively 
evolving galaxies, another possibility is that 
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they might be starbursting systems reddened 
by dust. Hence, an alternative explanation 
is that the quasars live in dense environments 
dominated by galaxies with enhanced star for- 
mation. In this case, the excess EROs may be 
dusty starburst galaxies, possibly the progen- 
itors of massive ellipticals in a forming clus- 
ter. There are examples of fields around high- 
redshift radio-loud AGN with excesses of star- 
forming galaxies such as Ly-alpha emitters 
(Kurk et al. 2000; Pentericci et al. 2000; Ven- 
emans et al. 2002). Also, Ivison et al. (2000) 
find a concentration of submm sources in the 
field of the z = 3.8 radio galaxy 4C 41.17 
raising the possibility that the EROs in this 
field are dusty starburst systems rather than 
passively evolving galaxies. This possibility is 
also discussed by Small et al. (2003) for a lu- 
minous submm source in the field of the radio 
galaxy 53W002 at z = 2.4. 

One might expect a correlation between the 
radio properties of the quasars and their ERO 
excesses if there is e.g. a relation between ra- 
dio power and black hole mass (Wold et al. 
2000; Laor 2000; Lacy et al. 2001). We have 
checked this for our sample, but do not find 
that the radio power correlates with the ERO 
excess. 

5.2. EROs in foreground lensing struc- 
tures 

Little is known yet of the luminosity func- 
tion of EROs. In the spectroscopic survey by 
Cimatti et al. (2002), which is 67 % complete 
at Ks < 19.2 mag, the R - Ks > 5.0 mag 
EROs are found to lie at z = 0.7 — 1.5. The 
redshift distribution derived by Miyazaki et 
al. (2002) based on photometric redshifts is 
in broad agreement with Cimatti et al's, but 
has a tail toward higher redshifts. If we as- 
sume that our ERO sample has a redshift dis- 
tribution similar to that found by Cimatti et 
al. and Miyazaki et al., we have to explain 
why the ERO surface density in the quasar 



fields is so significantly different from that in 
the large random-field surveys. We are sam- 
pling the ERO distribution with several small, 
but widely separated fields, and should there- 
fore be minimally affected by cosmic variance. 
Instead, the opposite seems to be the case; 
the quasar fields sample the densest regions in 
the ERO distribution. A large number of the 
quasar fields contain clumps of EROs, similar 
to what is seen only in the densest areas in the 
large random-field surveys. If the EROs have 
a peak redshift of z ^ 1.0, they must there- 
fore be clusters of galaxies, or overdensities in 
the ERO distribution, lying along the lines of 
sight to the quasars. 

A biasing of quasar lines of sight with re- 
spect to z !=a 1 overdensities in the ERO distri- 
bution can be explained by gravitational lens- 
ing. If we assume that the overdensities in the 
quasar fields do not extend much beyond the 
edges of the fields their typical scale is ~ 1 ar- 
cmin, corresponding to 2-3 Einstein radii for 
galaxy clusters at z ~ 1. The most efficient 
lens redshift for a quasar at z = 2 is 0.5-0.6, 
but a magnification can still occur if the lens 
lies at z ~ 0.8 - 1.0. 

The quasars in our sample are very lumi- 
nous, eight of them have My < —27.2 mag 
{Ho = 50 kms~^Mpc~^), and hence classify 
as highly luminous quasars. They therefore 
fulfill the three criteria for maximizing the 
probability of a magnification bias, i.e. they 
are distant, bright and have a steep number- 
magnitude slope (Claeskens Sz Surdej 2002, 
and references therein). If a magnification 
bias is the explanation for our results, we 
might expect that the fields around the most 
luminous quasars contain the highest overden- 
sities of EROs, but our data do not show any 
significant correlations between quasar lumi- 
nosity and ERO surface density to support 
this. 

Furthermore, if the majority of the EROs 
are associated with overdensities along the 
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line of sight to the quasars, we might ex- 
pect the spectra of the quasars in the rich- 
est fields to have strong (rest-frame equiv- 
alent width > 0.7 A) metal absorption line 
systems (Liu et al. 2000). We have found ab- 
sorption line data for seven of the 13 quasars 
in our sample (Junkkarinen, Hewitt, & Bur- 
bidge 1991, 1992). Two of the four richest 
fields have such data; Q2150+053 show no 
strong metal absorption lines and Q2338+042 
have CIV and Fell at Zabs ^1.80 (Barthel, 
Tytler, Sz Thomson 1990). Other quasars 
with metal absorption line systems in the fore- 
ground is Q0830-M15 with Fell at Zabs = 0.92 
and Mgll at Zabs = 0.80 (Sargent, Steidel, & 
Boksenberg 1989), Q0941-F261 with Mgll sys- 
tems at Zabs = 0.71 and Zabs = 1-09, and 
Q1354+258 with Mgll at Zabs = 0.86 and 
some weaker Fell and Mgll systems at Zabs = 
1.42 (Barthel, Tytler, & Thomson 1990), but 
these quasars lie in fields with relatively few 
EROs compared to the richer fields. The 
data for Q0758+120 and Q1629-f680 show no 
strong metal absorption line systems (Barthel, 
Tytler, t Thomson 1990). There is thus no 
clear trend of increasing number of absorp- 
tion line systems with ERO richness. Also, 
the quasars in our sample were not selected 
on the basis of having strong or numerous ab- 
sorption line systems. 

6. Conclusions 

On the average, the quasar fields stud- 
ied here have a cumulative surface density of 
EROs at Kg < 18, 18.5 and 19 mag which 
is two-three times larger than what is found 
in the general field. The excess is significant 
at the 3(7 level even after taking into ac- 
count that the ERO distribution is strongly 
clustered. There is also an excess of EROs 
at fainter magnitudes, but at Kg « 20 mag 
our counts become strongly affected by incom- 
pleteness and a secure quantitative estimate 
of the excess for Kg > 20 mag is difficult. 



However, by making a correction for the in- 
completeness, we find that the ERO excess at 
Ks < 19.5 and Kg < 20 mag may also be 
significant at the > 3a level. 

The quasar fields are small, typically 
2.1x2.1 arcmin^ (limited by the near- infrared 
field of view), but widely separated, and 
therefore sample the ERO population with 
minimal influence from cosmic variance. The 
fact that the excess is statistically signifi- 
cant therefore suggests that the quasar fields 
probe the densest regions in the ERO distri- 
bution, implying that the fields around lumi- 
nous, high-redshift quasars are on the average, 
overdense in EROs. There are two possible 
explanations for the observed overdensities. 
(1) The majority of the EROs may be as- 
sociated with concentrations of galaxies at 
the quasar redshifts, or (2) the EROs lie in 
the foreground, along the line of sight to the 
quasars. 

Based on both spectroscopic observations 
of z > 2 red galaxies (van Dokkum et al. 
2003) and spectral synthesis models (Fioc & 
Rocca-Volmerange 1997), we argue that a red- 
shift of z !=a 2 for the EROs is plausible. 
The Kg ~ 19 mag ERO excess may there- 
fore be physically associated with the quasars, 
in which case their magnitudes and colors are 
consistent with > L* passively evolving galax- 
ies. 

The EROs at i^s > 19 - 19.5 mag are 
more clustered in the immediate 20 arcsec 
region around the quasars, hence the radial 
distribution supports the explanation that 
the fainter EROs are physically associated 
with the quasars. Furthermore, the fraction 
of EROs among the JCg-selected galaxies in 
the quasar fields is unusually high, typically 
20-25 %, in contrast to < 10 % which is 
found in random-field surveys (Daddi et al. 
2000a; Scodeggio & Silva 2000; Roche et al. 
2002). This is consistent with the ERO ex- 
cess being associated with dense concentra- 
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tions of red, possibly evolved, galaxies. The 
high ERO fractions may therefore indicate 
that a morphology-density relation (Dressier 
et al. 1997) exists in the quasar fields. The 
candidate ERO clusters in the quasar fields 
may therefore represent the oldest and most 
evolved systems at z > 2, and will be crucial 
for testing models of structure and galaxy for- 
mation. If the ERO excess is dominated by 
dusty starforming galaxies instead of passively 
evolving galaxies, the EROs might be the pro- 
genitors of massive ellipticals in a forming 
cluster associated with the quasars. 

The second explanation, that the ERO ex- 
cess lies in the foreground, is equally interest- 
ing because it implies that the lines of sight 
toward high-redshift, luminous quasars are bi- 
ased toward overdensities in the ERO distri- 
bution. In this case, the ERO excess is most 
likely associated with clusters of galaxies at 
2; pa 1 giving rise to a magnification bias of the 
quasars. In this interpretation, we would be 
seeing candidate clusters with passively evolv- 
ing galaxies at z ?a 1, a redshift where still 
only less than a handful of clusters are known. 
However, we see no trend of quasar brightness 
with ERO overdcnsity as might be expected in 
this interpretation. Also, there is no associa- 
tion of the numbers of absorption line systems 
with ERO richness for seven of the 13 quasars 
with absorption line data. 

With the current data, it is not possible to 
tell whether the ERO excess is part of a small- 
scale overdensity centered on, or lying close to, 
the quasars, or if it is part of a larger structure 
that extends over a wider area. To address 
this issue, we have begun an imaging program 
using J- and Kg-filters at the Palomar 200- 
inch telescope over much larger fields. With 
these data we will be able to trace the ERO 
population around the z '^1 quasars on scales 
up to 4-5 Mpc. 

We are grateful to Dave Thompson, Lin 
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tute of Technology, funded by the National 
Aeronautics and Space Administration and 
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search has also made use of the NASA/IPAC 
Extragalactic Database (NED) which is oper- 
ated by the Jet Propulsion Laboratory, Cali- 
fornia Institute of Technology, under contract 
with the National Aeronautics and Space Ad- 
ministration. 
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